Galerkin and Runge–Kutta methods: unified formulation, a posteriori error estimates and nodal superconvergence

Abstract. We unify the formulation and analysis of Galerkin and Runge–Kutta methods for the time discretization of parabolic equations. This, together with the concept of reconstruction of the approximate solutions, allows us to establish a posteriori superconvergence estimates for the error at the nodes for all methods. 1

Antenatal Sonographic Diagnosis of Pharyngeal Teratoma: Our Experience of a Rare Case with Review of the Literature

Background. Teratomas are the most common tumors. They are usually localized in the sacrococcygeal area, while the pharyngeal localization is very rare. The number of cases of stomatopharyngeal teratomas detected prenatally via sonography is very small. Case Report. We present the case of a 24-year-old primipara at 18 weeks' gestation, that at the routine ultrasound scan, the fetus was found with an echogenic mass, filling the stomatopharyngeal cavity and protruding from the mouth. Other abnormalities were not found. Termination of pregnancy was achieved using misoprostol. A female stillborn fetus with a weight of 250 g and length of 25.5 cm was delivered. The postmortem and pathologic examination confirmed the diagnosis. Conclusion. Pharyngeal teratomas can be diagnosed with the use of ultrasounds in utero facilitating parents' counseling in early time

A system of ODEs for a Perturbation of a Minimal Mass Soliton

We study soliton solutions to a nonlinear Schrodinger equation with a saturated nonlinearity. Such nonlinearities are known to possess minimal mass soliton solutions. We consider a small perturbation of a minimal mass soliton, and identify a system of ODEs similar to those from Comech and Pelinovsky (2003), which model the behavior of the perturbation for short times. We then provide numerical evidence that under this system of ODEs there are two possible dynamical outcomes, which is in accord with the conclusions of Pelinovsky, Afanasjev, and Kivshar (1996). For initial data which supports a soliton structure, a generic initial perturbation oscillates around the stable family of solitons. For initial data which is expected to disperse, the finite dimensional dynamics follow the unstable portion of the soliton curve.Comment: Minor edit

hp-adaptive Galerkin Time Stepping Methods for Nonlinear Initial Value Problems

This work is concerned with the derivation of an a posteriori error estimator for Galerkin approximations to nonlinear initial value problems with an emphasis on finite-time existence in the context of blow-up. The structure of the derived estimator leads naturally to the development of both h and hp versions of an adaptive algorithm designed to approximate the blow-up time. The adaptive algorithms are then applied in a series of numerical experiments, and the rate of convergence to the blow-up time is investigated

Finite Element Convergence for the Joule Heating Problem with Mixed Boundary Conditions

We prove strong convergence of conforming finite element approximations to the stationary Joule heating problem with mixed boundary conditions on Lipschitz domains in three spatial dimensions. We show optimal global regularity estimates on creased domains and prove a priori and a posteriori bounds for shape regular meshes.Comment: Keywords: Joule heating problem, thermistors, a posteriori error analysis, a priori error analysis, finite element metho

Second order averaging for the nonlinear Schroedinger equation with strongly anisotropic potential

International audienceWe consider the three dimensional Gross-Pitaevskii equation (GPE) describing a Bose-Einstein Condensate (BEC) which is highly confi ned in vertical z direction. The highly confi ned potential induces high oscillations in time. If the confi nement in the z direction is a harmonic trap (which is widely used in physical experiments), the very special structure of the spectrum of the confi nement operator will imply that the oscillations are periodic in time. Based on this observation, it can be proved that the GPE can be averaged out with an error of order of epsilon, which is the typical period of the oscillations. In this article, we construct a more accurate averaged model, which approximates the GPE up to errors of order epsilon squared. Then, expansions of this model over the eigenfunctions (modes) of the vertical Hamiltonian Hz are given in convenience of numerical application. Effi cient numerical methods are constructed for solving the GPE with cylindrical symmetry in 3D and the approximation model with radial symmetry in 2D, and numerical results are presented for various kinds of initial data

A posteriori error estimates for fully discrete schemes for the time dependent Stokes problem

The final publication is available at Springer via http://dx.doi.org/10.1007/s10092-018-0259-2This work is devoted to a posteriori error analysis of fully discrete finite element approximations to the time dependent Stokes system. The space discretization is based on popular stable spaces, including Crouzeix–Raviart and Taylor–Hood finite element methods. Implicit Euler is applied for the time discretization. The finite element spaces are allowed to change with time steps and the projection steps include alternatives that is hoped to cope with possible numerical artifices and the loss of the discrete incompressibility of the schemes. The final estimates are of optimal order in L∞(L2) for the velocity error